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Abstract 

We investigate how well the shape parameterizations are applicable for studying the 
giant dipole resonance (GDR) in nuclei, in the low temperature and/or high spin 
. regime. The shape fluctuations due to thermal effects in the GDR observables are 

calculated using the actual free energies evaluated at fixed spin and temperature. 
The results obtained are compared with Landau theory calculations done by param- 
eterizing the free energy. We exemplify that the Landau theory could be inadequate 
where shell effects are dominating. This discrepancy at low temperatures and high 
spins are well reflected in GDR observables and hence insists on exact calculations 
in such cases. 
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In recent years considerable interest has been shown [1,2,3,4] to study the 
structural transitions as a function of both angular momentum and temper- 
ature in highly excited nuclei. The Giant Dipole Resonance (GDR) studies 
have been proved to be a powerful tool to study such hot and rotating nuclei 
[5] and recently the domain of GDR spreads rapidly over different areas of 
theoretical and experimental interest [6,7,8,9]. The GDR observations provide 
us information about the geometry as well as the dynamics of nuclei even 
at extreme limits of temperature (T), spin (/) and isospin (r). In the past 
most of the GDR measurements in hot nuclei were made at moderate and 
high T [5,6,7,8,9,10,11,12]. Few experiments have been carried out recently to 
study the GDR states at low temperatures [1,2,3]. Hence the theories which 
were successful in the high T regime should be now scrutinized with the low 
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temperature observations as well as with theories incorporating properly the 
microscopic effects (such as shell effects) which are dominant at low tem- 
peratures. In a macroscopic approach, while dealing with the thermal shape 
fluctuations, free energy parameterizations such as Landau theory [13,14,15] 
are usually employed to do timesaving calculations. In this work we survey 
the applicability of Landau theory by demanding consistency with exact cal- 
culations done without any parameter fitting. 

The theoretical approach we follow is of three fold with models for 1) shape 
calculations, 2) relating the shapes to GDR observables and 3) considering 
the shape fluctuations due to thermal effects. For shape calculations we follow 
the Nilsson-Strutinsky (NS) method extended to high spin and temperature 
[16]. The total free energy (Ftot) at fixed deformation is calculated using the 
expression 

F TO t = -Erldm + J~]5F ■ (1) 

p,n 

Expanding the rotating liquid-drop energy -Erldm and writing shell corrections 
in rotating frame [17] leads to 

£ldm + E 5F " + iU( 'tot + E 6I ) ■ ( 2 ) 

p,n ^ p,n 

The angular velocity u is tuned to obtain the desired spin given by 

/ TOT = <3 rigUJ + SI . (3) 

The liquid-drop energy (E LDM ) is calculated by summing up the Coulomb and 
surface energies [18] corresponding to a triaxially deformed shape defined by 
the deformation parameters /3 and 7. The rigid-body moment of inertia (Qvig) 
is calculated with surface diffuseness correction [18]. The shell correction (5F U ) 
is the difference between the deformation energies evaluated with a discrete 
single-particle spectrum and by smoothing (averaging) that spectrum and is 
given by the relation 

Ee^-T$> )-(E^- T E*0 > ( 4 ) 

i=l J \ i i J 

where the discrete quantities 

1 

1 + exp [- jl t-) 

Si = — [rii \nrii — (1 — n,j) ln(l — n*)], the Strutinsky smeared quantities n« = 
f^°ocf( x ) n i( x ) dx, Si = Jf^, f(x) Si(x) dx and f(x) is the smearing function 
[16]. Similarly the shell correction corresponding to the spin is given by 

00 00 

SI = I - I = E m » n i ~ E m ^ • ( 6 ) 
i=l i=l 
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The single-particle energies (ef ) and spin projections (m,) are obtained by 
diagonalizing the triaxial Nilsson Hamiltonian in cylindrical representation 
upto first twelve major shells. 

In a macroscopic approach, the GDR observables are related to the nuclear 
shapes. This is realized using a model [19] comprising an anisotropic harmonic 
oscillator potential with separable dipole-dipole interaction. In this formalism 
the GDR frequencies in laboratory frame are obtained as 

u t = (l + r,)^u; g , (7) 
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^3 t n = I (1 + u)-^ + & - 2 [(1 + ^) 2 K - ^) 2 

+8fi 2 (l +?7 )(^ + ^}%Q, (9) 

where f2 is the cranking frequency, w x , lo v , lo z are the oscillator frequencies 
derived from the deformation of the nucleus and 77 is a parameter that charac- 
terizes the isovector component of the neutron and proton average field. The 
GDR cross sections are constructed as a sum of Lorentzians given by 

a{E 1 )=T\ ^—o , (10) 

where Lorentz parameters E m , a m and T are the resonance energy, peak cross- 
section and full width at half maximum respectively. Here % represents the 
number of components of the GDR and is determined from the shape of the 
nucleus [19,20]. The energy dependence of the GDR width can be approxi- 
mated by [21] 

Ti « 0.02QE1 9 . (11) 
The peak cross section a m is given by 

°m = 6 °-^-f 0.86(1 + a) . (12) 

The parameter a which takes care of the sum rule is fixed at 0.3 for all the 
nuclei considered in this work. In most of the cases we normalize the peak 
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with the experimental data and hence the choice of a has less effect on the 
results. The other parameter 77 varies with nucleus so that the ground state 
GDR centroid energy is reproduced. The choice for 84 Zr is rj = 2.6, for 147 Eu, 
rj = 2.8, for 179 Au, r] = 3.2, and for 208 Pb, rj = 3.4. For calculating the GDR 
width, only the power law (11) is used in this work and no ground state width 
is assumed. 

When the nucleus is observed at finite excitation energy, the effective GDR 
cross-sections carry information on the relative time scales for shape rearrange- 
ments [22] which lead to shape fluctuations. In the case of hot and rotating 
nuclei, apart from thermal shape fluctuations, there can be fluctuations in the 
orientation of the nuclear symmetry axis with respect to the rotation axis. The 
general expression for the expectation value of an observable O incorporating 
both thermal and orientation fluctuations is given by [13,23] 

where Q = (0, 9, ip) are the Euler angles specifying the intrinsic orientation 
of the system, Cj -X • Cj — Sy^ cos 2 sin 2 9 + 3y</ sin 2 sin 2 9 + Qy z / cos 2 9 
is the moment of inertia about the rotation axis Co given in terms of the 
principal moments of inertia 9y x ', ^y'y 1 -, ^z>z>, and the volume element T>[a] = 
(3 4 1 sin 37 1 df3 g?7 sin 9d9d<p. 

The study of thermal fluctuations by numerical evaluation of Eq. (13) in gen- 
eral requires an exploration of five dimensional space spanned by the deforma- 
tion and orientation degrees of freedom, in which a large number of points are 
required in order to assure sufficient accuracy (especially at finite angular mo- 
mentum). Hence certain parameterizations were developed [14,24] to represent 
the free energy using functions that mimic the behaviour of the NS calcula- 
tion as closely as possible. One such parametrization is the Landau theory of 
phase transitions, developed by Alhassid and collaborators [13,14,22,23]. Here 
the free energy is expanded in terms of certain temperature dependent con- 
stants which are to be extracted by fitting with the free energy calculations 
at fixed temperatures from the NS method. Moreover, once the fits involving 
free energy and moment of inertia are made for the non-rotating case, the 
calculations can be extended to higher spins using the relation [13] 

F(T, J; (3, 7, A) = F(T, u = 0; A 7) + \ 7 ' K ■ (14) 

z <jj • J, ■ uo 

Hence this theory offers an economic parametrization to study the hot rotating 
nuclei. However the above expression carries the shell corrections evaluated at 
uj = all along to higher spins. This is not desirable as the single-particle levels 
swiftly change positions with increasing spin, resulting in a totally different 
shell structure. We have employed Landau theory in its extended form as given 
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in Refs. [14,16]. 



Recently [16,25] few calculations have been done by performing the thermal 
fluctuation calculations exactly by computing the integrations in Eq. (13) 
numerically with the free energies and the observables being calculated "ex- 
actly" at the integration (mesh) points. In this way the calculations can be 
done more accurately without using any parametrization and consequent fit- 
ting. In this work we have performed such calculations, however, neglecting 
the orientation fluctuations. This enables us to perform the integration in the 
deformation space only which at present is two dimensional having the defor- 
mation parameters f3 and 7. The role of orientation fluctuations is negligible 
while calculating the scalar observables [16,26] such as the GDR cross section 
and width. The derivation [26] of Eq. 13 employs the assumption of Eq. (14) 
and here we discuss the consequences as the free energy now takes the form 
of Eq. (1). The partition function at fixed angular momentum I is obtained 
as [26] 



(15) 



The exponent in the above integral can be taken as —F^qt/T in our case. 
Letting Co - X • Cj = Qyy = Stot = ^rig + and following the definitions 
given in Ref. [17], we have 

d J e -F TOT /Tj = c- TQT u e -F TOT /T duJ _ 

In the limit g^tot / doo — > 0, the integration over 00 in Eq. (15) is analytically 
solvable leading to an expression having exactly same form as Eq. (13) with 
Co ■ X ■ Co replaced by 3tot- The shell corrections to moment of inertia are 
significant for spherical shapes which are however suppressed by the term /3 4 
in V[a] and for (3 > 0.1 we have Stot ~ ^rig- At high spins also the limiting 
condition is very much valid as the well deformed shapes are more favoured. 
It has already been shown [26] that this transformation from frequency to 
angular momentum leads to a pre-f actor to the volume element of the integral. 
However we show later that this factor has not much role to play in the 
practical calculations. While performing fluctuation calculations in this way, 
the free energy at any given spin is obtained by tuning the cranking frequency 
to get the desired spin. In such case the projection is not necessary and has 
been neglected in some similar calculations 



Now we compare our calculated results obtained by using the extended Lan- 
dau theory and the exact method. The calculations are performed with 1) 
the liquid-drop model (LDM) free energies and Landau theory, 2) NS free 
energies and Landau theory and 3) NS free energies with exact treatment of 
fluctuations. 
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Fig. 1. GDR width for the nuclei 208 Pb and 179 Au. The results obtained using liquid 
drop model (dotted line), Landau theory (dashed line) and the exact calculations 
(solid line) are compared. Left: Experimental data represented by solid squares are 
taken from ref. [11] and the revised data [27] are represented by solid triangles. 
Right: The different curves carry same meaning as in the left panel. At T = 
MeV, the widths calculated assuming oblate (fa = —0.22), prolate (/% = 0.24) 
and spherical shapes are denoted by the open symbols down-triangle, up-triangle 
and circle respectively. The experimental value [2] at T = MeV is the intrinsic 
(spherical) width Tq and at T = 0.7 MeV is the total width having To = 5 MeV 
and 02 = 0.1. Both these values are represented by solid circles. 
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Fig. 2. Potential energy surfaces for 179 Au at T = 0.0, 0.7 and 1.0 MeV, calculated 
using the Nilsson-Strutinsky method. The contour line spacing in 0.5 MeV, the most 
probable shape is marked by a solid circle and the first two minima are shaded. The 
barrier between the co-existing shapes diminish with increase in temperature and 
the nucleus becomes 7-soft at T ~ 0.7 MeV. 

In Fig. 1 we present the calculated GDR widths of 208 Pb and 179 Au at u = 
along with experimental results. In the case of 208 Pb, strong shell corrections 
for spherical shape results in large difference in the deformation energies be- 
tween spherical and deformed configurations. This leads to attenuation of 
thermal fluctuations at lower temperature and hence the obtained widths are 
much lower when compared to liquid-drop model results. The magnitude of 
this attenuation comes out to be different in methods 2 and 3. More discus- 
sions and comparison with other reported results can be found in Ref. [16]. 
Recently [2] the thermal fluctuation calculations of 208 Pb were extrapolated 
to interpret low T measurement of in 179 Au. This extrapolation was as- 

sumed to give the lower limit for the T GDR . Our calculations as shown in Fig. 
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1 do not favour any such interpretation. In fact the situation in 179 Au is totally 
different from that of 208 Pb as now deformed shapes are more favoured. The 
ground state deformation comes out in our calculations as f3 2 — —0.22 which is 
slightly larger than f3 2 — —0.17 obtained by relativistic mean field theory (with 
NL3 parameter set) and Hartree-Fock-Bogoliubov calculations [28]. T GDR cal- 
culated for this deformation at T = MeV lie very close to those obtained 
with exact thermal fluctuations at T = 0.1 MeV. Moreover coexisting oblate 
and prolate shapes also prevail as we observe in the potential energy surfaces 
(PES) given in Fig. 2. We infer from the PES that the deformed shapes only 
are favoured even at T = 0.7 MeV where we see 7-softness. As the T GDR of 
deformed shapes are always higher than spherical ones, the calculation with- 
out shell effects (LDM), which favour spherical shapes, give the lower limit 
for Tgdr i n deformed nuclei. Hence the inclusion of shell corrections will only 
increase T GDR and could not account for the suppression of T GDR in 179 Au. 

A similar situation has been observed in 120 Sn [1] also where the thermal 
fluctuation model, even with shell effects included, could not explain the sup- 
pression of T GDR at low T. Our exact calculations (see Fig. 4 of Ref. [16]) also 
could not explain this anomaly. It has been speculated [29] that this suppres- 
sion could be a general feature of all nuclei, independent of shell effects. The 
possible cause of this could be the pairing correlations which can be significant 
at low T. The inclusion of fluctuations in the pairing field [30] or fully micro- 
scopic calculations [31] considering the temperature dependence of spreading 
width could explain the low width of 120 Sn at low T. The same could be the 
case of 179 Au also. In the present work, we wish to emphasis on the problems 
in using the shape parameterizations apart from these effects. 

Also we infer from the calculations of 208 Pb and 179 Au, if the shell effects are 
strong leading to a crisp or multiple minima, the Landau theory could not 
account for it. This discrepancy can be partially ascribed to the parameteri- 
zation itself and the fitting procedure involved. It was already suggested [14] 
that at low temperatures the least square fit for Landau constants could be 
erroneous and techniques like uniform mapping should be adopted. The new 
parameterization [24] may solve some of these issues despite having some in- 
consistencies at higher temperatures. However in both the methods involving 
shape parameterization, the high spin calculations are done using Eq. (14). 
Due to the complexity of shell structure at high spins, the difference in shell 
effects with changing spin may be crucial especially at lower temperatures. 

In our previous work [16] we investigated the spin dependence of T GDR in 
120 Sn at T = 1.8 MeV (See Fig. 5 of Ref. [16]) and found the Landau theory 
to perform well even at spins upto 70h. This is understandable as there is no 
spin dependence of shell effects in this case and the temperature is sufficiently 
high where the Landau theory works well. Our results for 179 Au at T = 0.7 
MeV and at different spins are shown in Fig. 3. Here also the spin dependence 
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Fig. 3. Spin dependance of GDR width in Au and Eu. The results obtained 
using liquid drop model (dash-dotted line), Landau theory (dashed line) and the 
exact calculations (solid line) are compared. Right: The solid circle, solid square 
and solid triangle correspond to experimental data [12] at beam energies 170, 165 
and 160 MeV respectively. These energies correspond to temperatures from 1.2 to 
1.4 MeV [12]. 

of shell effects are found to be insignificant. We extended our study to 147 Eu 
for which experimental data is available for T ~ 1.3 MeV. Our results are 
shown in Fig. 3 where we can see that at T = 1.3 MeV there is not much 
difference between the results of Landau theory and exact calculations. Also 
the widths are very much similar to those obtained using LDM as the proton 
and neutron shell corrections are weak and they act against themselves. This 
trend continues even at spins up to 60/1. However at T = 0.5 MeV, the situation 
is drastically different as the spin-driven shell effects play their role. At uo = 0, 
the shell correction is of the order of 2 MeV and hence the three methods 
give different results. The shell corrections decrease with the increasing spin, 
at 40H the equilibrium shape acquires deformation 7) = (0.3,-180°) and 
at 60h, it is (0.3,-130°). The effect of these sharp changes survive thermal 
fluctuations in the exact calculations and is averaged out in the Landau theory 
as well as the LDM calculations. It has to be noted that the Landau theory 
results carry the shell corrections calculated at uo = all along. From the above 
arguments it is clear that one cannot substantiae the success of Landau theory 
in the cases of 120 Sn and 179 Au. For certain nuclei like 147 Eu the spin-driven 
shell effects can be crucial at low T. 

In Fig. 4 we show for the hot rotating 84 Zr nucleus, the results of our GDR cross 
section calculations. It has been observed that in neutron-deficient Zr isotopes, 
spin-driven shell effects are stronger leading to a sharp shape transition at 
lower temperatures [16]. It is evident from the Fig. 4 that, in certain nuclei 
at higher spins the spin-driven shell effects may play vital role even at T ~ 1 
MeV. Also we show in Fig. 2 the impact of including the pre-factor ^tot 
in the thermal fluctuation integral. This is the only case presented in this 
work where we could notice a slight effect of the factor. One can observe that 
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Fig. 4. GDR cross sections for the nuclei 84 Zr at T = 1.0 MeV and / = 60h with 
shape fluctuations using Landau theory (dashed line) , exact method with the factor 
^tot (Exact 1, solid black line) and exact method without that factor (Exact 2, 
solid grey line). 

even at extreme limits of spin the pre-factor has practically no effect as the 
dominating role is played by the exponential term exp(— F^qt/T) which has 
exact temperature and spin dependence. 

To summarize, in this work the thermal fluctuations are dealt in an exact way 
without any parameter fitting. Comparison of our present approach with the 
thermal fluctuation model comprising Landau theory suggests that the shape 
parameterizations could be insufficient for GDR calculations in the presence 
of strong shell effects. The discrepancies are shown in certain nuclei at two 
situations namely at very low temperatures (T < 1 MeV) and at moderate 
temperatures and high spin (T ~ 1 MeV, I > 30h). The latter case is ascribed 
to spin driven shell effects, which the existing shape parametrizations do not 
account for. The present study necessitates exact treatment of fluctuations in 
these cases where there is experimental and theoretical focus recently. 
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